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K* ' Abstract 

CS ; Recent T2K, MINOS and Double CHOOZ results, together with global fits of 

mixing parameters, indicate a sizeable reactor angle ^13 ~ 8° which, if confirmed, 
would rule out tri-bimaximal (TB) lepton mixing. Recently two of us studied 

f^ , the vacuum alignment of the Altarelli-Feruglio A^ family symmetry model includ- 

CN ' ing additional flavons in the 1' and 1" representations, leading to so-called "tri- 

maximal" neutrino mixing and allowing a potentially large reactor angle. Here 
we show how such a model may arise from a Supersymmetric (SUSY) Grand 

^ \ Unified Theory (GUT) based on SU{5), leading to sum rule bounds \s\ < -^, 

H I |a| < ^(r + -^)| cos S\, where s, a, r parameterise the solar, atmospheric and reac- 

tor angle deviations from their TB mixing values, 6 is the CP violating oscillation 
phase, and 6c is the Cabibbo angle. 
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1 Introduction 

Recently T2K have published evidence for a large non-zero reactor angle [1]. Combining 
this with data from MINOS [2] and other experiments in a global fit yields [3] 613 = 
6.5° ± 1.5°, assuming a normal neutrino mass hierarchy. Here the errors indicate the 
la range, although the statistical significance of a non-zero reactor angle is about 3cr. 
Other global fits indicate a larger reactor angle 9is = 9.1° ± 1.3° |l]|j The first results 
from Double CHOOZ also show indications of a non-zero reactor angle [6]. If confirmed, 
a sizable reactor angle 6*13 ~ 8°, consistent with both global fits, would rule out the 
hypothesis of exact tri-bimaximal (TB) mixing [7]. 

In the framework of Supersymmetric (SUSY) Grand Unified Theories (GUTs) of 
Flavour [S] (i.e. with a Family Symmetry [5] implemented) it is already known that TB 
mixing cannot be exact. The reason is that, even if TB mixing is realised exactly in the 
neutrino sector, observable lepton mixing is subject to charged lepton (CL) corrections 
(due to the fact that L^pmns = UeUl) and renormalisation group (RG) corrections, not to 
mention other corrections due to canonical normalisation (CN) (for a unified discussion 
of all three corrections see e.g. [10] and references therein). Therefore, in the framework 
of SUSY GUTs of Flavour, the question of whether TB mixing may be maintained in 
the neutrino sector is a quantitative one: can the above CL, RG and CN corrections be 
sufficiently large to account for the observed reactor angle? The answer is yes in some 
cases (see e.g. [H]), but no in many other cases. For example, in models based on the 
Georgi-Jarlskog mechanism [12], where the CL corrections are less than or about 3°, 
and where the RG and CN corrections are less than or about 1° (which is the case for 
hierarchical neutrinos), it would be difficult to account for a reactor angle 6'i3 ~ 8°. For 
this reason, there is a good motivation to consider other patterns of neutrino mixing 
beyond TB mixing, and many alternative proposals [13] have indeed been put forward to 
account for a non-zero Oi^. On the other hand, since the solar and atmospheric mixing 
angles remain consistent with TB mixing, there is also a good motivation to maintain 
these successful predictions of TB mixing. 

In going beyond TB mixing, it is useful to consider a general parameterisation of the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix in terms of deviations from 
TB mixing 



^PMNS 



^^{l + s-a + re'') ^^{l - \s - a ~ \re^') ^(1 + a 
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\ ^^[l + s + a-re^') --i=(l - 1, + a + Ire*^) ^(1 - a) 

where the deviation parameters s, a, r are defined as [H], 

1 1 r 

sin6'i2 = ^(1 + s) , sin6'23 = ^(1 + a) , sin 6^13 = —= . (2) 



^In both cases, the quoted ranges are calculated using the new reactor anti- neutrino fluxes [5]. 



For example, the global fit of the conventional mixing angles in [3] can be translated 
into the la ranges: 

0.12 < r < 0.20, -0.06 < s < -0.006, -0.06 < a < 0.08. (3) 

In a SUSY GUT of Flavour, the Family Symmetry is responsible for determining the 
neutrino mixing pattern, which then gets corrected by CL, RG and CN contributions 
to yield the observed lepton mixing angles. The question is what is the underlying 
neutrino mixing pattern? If we want to go beyond TB neutrino mixing, there are many 
possibilities. One simple scheme is the trimaximal (TM) mixing pattern |15j : 

/ ^cos^ ^ ^sin^e^^ \ 

UtI = P' -^cos^-^sin^e-" ^ ^ cos^ - ^ sin^e^H P , (4) 
\^__i=cos^+^sin^e-*^ 73 -^ cos^ - ^ sin^e^^y 



where -Tgsin'i? = sin6'5i'3, P' is a diagonal phase matrix required to put ?7pmns = U^U^^Ia 

into the PDG convention [Id], and P = diag(l, e*^", e*^") contains the usual Majorana 
phases. In particular TM mixing approximately predicts TB neutrino mixing for the 
solar neutrino mixing angle 612 ~ 35° as the correction due to a non-zero but relatively 
small reactor angle is of second order. However we emphasise again that, in a SUSY 
GUT of Flavour, TM mixing refers to the neutrino mixing angles only, and the physical 
lepton mixing angles will involve additional CL, RG and CN corrections. Nevertheless, 
TM neutrino mixing could provide a better starting point than TB neutrino mixing, 
if the reactor angle proves to be 613 ~ 8°, and this provides the motivation for the 
approach followed in this paper. 

Recently, an A4 model of TM neutrino mixing was discussed in [17]. In the original 
A4 models of TB mixing Higgs fields [12] or flavon fields [19J transforming under A4 as 
3 and 1 but not 1' or 1" were used to break the family symmetry and to lead to TB 
mixing. However there is no good reason not to include flavons transforming as 1' or 1", 
and once included they will lead to deviations from TB mixing [20], in particular it was 
noted that they lead to TM mixing [21] . In [17] the vacuum alignment of the Altarelli- 
Feruglio A4 family symmetry model [19], including additional flavons in the 1' and/or 
1" representations, was studied and it was shown that it leads to TM neutrino mixing. 
In this paper we shall show how such a model may arise from a SUSY GUT based on 
SU{5), leading to the sum rule bounds \s\ < ^ and \a\ < |(r + ^)| cos5|, up to RG 
and CN corrections, where r, s, a are the above tri-bimaximal deviation parameters, S 
is the CP violating oscillation phase, and 6c is the Cabibbo angle. However we shall 
not be interested in the details of the GUT breaking Higgs potential, which is model 
dependent and may not even be describable by renormalisable 4d field theory, as is the 
case where GUT breaking is due to orbifold constructions [22] . Instead we shall take the 
GUT breaking as "given" and formulate the theory just below the GUT scale, although 
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Table 1: Matter and Higgs chiral superfields in tlie model. 

we shall use SU{5) notation for convenience, and to emphasise the phenomenological 
predictions which follow from GUTs. 

The rest of the paper is organised as follows. In Section [2] we introduce the model, 
presenting field content, charges, flavon alignments and leading order (LO) superpoten- 
tial terms. Section E] then presents the mass matrices and mixing angles for neutrinos, 
quarks and charged leptons arising from the LO superpotential. The effect of the non- 
trivial charged lepton corrections (due to the grand unified setup) on the physical lepton 
mixing angles is discussed in Section HI We conclude in Section [51 The discussion of the 
vacuum alignment and the next to leading order (NLO) terms is presented in Appendices 
RJandlBl respectively. 

2 The model 



The transformation properties of the SU{5) matter and Higgs multiplets are shown in 
Table [TJ A^ denotes the right-handed neutrino fields, F the 5 of SU{5), containing 
the lepton doublet and the right-handed down- type quark, and T labels the 10 which 
includes the quark doublet as well as the right-handed up-type quark and charged lepton. 
A^ and F furnish the triplet representation of A4, thus unifying the three families of 
leptons, while the three families of the Tj transform in the three distinct one-dimensional 
representations of A^. In the Higgs sector, we have introduced the ii/45 in order to 
implement the Georgi-Jarlskog mechanism [T2]Jn 

The full set of flavon fields is shown in Table [2j The fields ips and ^* are responsible 
for the flavour structure of the neutrino sector, while the fiavons ifT and 6^ control the 
quark and charged lepton sector. The vacuum structure is obtained via the standard 



^The standard MSSM /x-term jiHuHd (where Hu is the SU{2)l doublet of Hs; and Hd is a linear 
combination of the SU{2)l doublets in H-g and -H45) is forbidden by the A4, U{1), Z3 and Z5 symmetries 
as well as U{1)r, allowing for a natural solution to the /^-problem of the MSSM using a GUT singlet 
from the hidden sector of Supergravity theories [23] . 
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Table 2: Flavon chiral superfields in the model. 

F-term alignment mechanism [19] where the F-terms of so-called driving fields (denoted 
by a superscript 0) are set to zero, thus giving rise to constraints which in turn fix the 
fiavon alignments. As shown in Appendix [Xj one obtains the following triplet fiavon 
alignmentsjj 

{Vt) oc IO , (^5) oc (ij. (5) 

The f/(l)_R represents an /^-symmetry whose Z2 subgroup gives rise to the standard 
i?-parity which forbids unwanted operators contributing to proton decay and keeps the 
lightest SUSY particle a good candidate for cold dark matter. The U{1) and the three 
Zjv shaping symmetries constrain the structure of the Yukawa matrices in the quark 
and charged lepton sectors. Specifically, the Z5 prevents the neutrino flavons {ips and 
^^) from appearing in the quark and charged lepton Yukawa couplings. 

In the neutrino sector, the ^4 family symmetry is broken by the fiavon fields ips 
and .^*, thereby leading to a TM mixing pattern as observed in [17]. In the quark and 
charged lepton sector the A^ symmetry is broken differently by virtue of the fiavon 
fields ipT and 6**. Due to the SU{b) structure, the form of the charged lepton and down 
quark Yukawa matrices is intimately related, leading to a non-trivial left-handed charged 
lepton mixing which combines with the TM structure of the neutrino mixing to give the 
physical PMNS mixing. 



2.1 Allowed terms 

The neutrino sector is composed of Dirac and Majorana mass terms which take the 
leading order form in the superpotential. 



W, = yFNH, + {y^^s + 2/2^ + Z/a^ + 2/2^") NN 



(6) 



■^The auxiliary fiavon field a is introduced for the purpose of achieving the alignment of the U{1) 
charged fiavon field (pr- 



with y, 2/1, j/25 2/25 y'i being dimensionless couplings. 

The leading order superpotential terms of the down quark and charged lepton sector 
are given as follows 

+ ^^^ (F^t) H^n + (^ (FcpT)' + ^ (F^t)"^ H^T^ (7) 



iFipT) + ^ {{FifT)") j i/5^3 + {^{FifT)'^ H^n 



where A^; is the relevant messenger mass. Note that for some entries of the down quark 
Yukawa matrix, there are several different operators of the same order; here we simply 
choose an example for illustrative purposes. The flavons 6'* play a role similar to a 
Froggatt-Nielsen field [23]. 

Finally the leading order up quark sector Yukawa superpotential terms take the form 

9^9'? (9^{9'f{9"f a9(9'?9'\, 

\ U^ U J ^g^ 

9^9' , 99' 9'(9")\ 

u u u 

It should be mentioned that the messenger mass in this sector, A„, may in principle be 
different from that in the down quark sector. The field 9' is introduced specifically to 
generate the T2T2 term to the required order. 

Examples of the many subleading higher order operators allowed by the symmetries 
of the model are listed in Appendix |Bl As their contribution to the mass matrices is 
negligible, they do not induce physically relevant modifications of the LO picture. 

3 Fermion mass matrices 

After spontaneous breakdown of the A4 family symmetry by the flavon VEVs, the 
superpotential terms of Eqs. ([S])-® predict mass matrices for the respective sectors. In 
the following, order one coefficients in the quark and charged lepton sectors are omitted 
(including flavon VEV normalisation factors) . Regarding the scale of the flavon VEVs 
we define 

'n = ^. (9) 

where ipi=ipT, 9^ or a. In order to get the hierarchical structure of the quark and charged 
lepton mass matrices we assume 

r]p = e^ and r^others = e, (10) 
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where the numerical values for e depend on the messenger scale of the relevant sector. 
We present LO operators for each entry in the mass matrices; NLO operators can be 
found in Appendix [Bl 

In the Higgs sector, it is not the if 5, H-^ or H-^ which get VEVs but their SU{2)l 
doublet components. These are the two MSSM doublets Hu (corresponding to H5) and 
Hd (corresponding to a hnear combination of H-^ and H^); they originate below the 
GUT scale and remain massless down to the electroweak scale. The non-MSSM states all 
acquire GUT scale masses, including the linear combination of H-g and Hjg orthogonal 
to Hd- Electroweak symmetry is broken after the light MSSM doublets Hu^ acquire 
VEVs Vu,d and they then generate the fermion masses. 

In the following all quark and charged lepton mass matrices are given in the L-R 
convention, i.e. the mass term for a field tp is given in the order i/jlMljiijjr. 

3.1 Neutrino sector 

Eq. (^ gives Dirac and Majorana mass matrices 
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(12) 



with a = yii^s), (3 = 2/2(0, Y = 2/2 (f) and 7" = 2/2 (D- As shown in [H], the 
standard type I seesaw formula then yields a light neutrino mass matrix of trimaximal 
structure, and hence a neutrino mixing matrix of the form as given in Eq. (jlj). 

3.2 Down quark and charged lepton sector 

In the down quark and charged lepton sector, the superpotential of Eq. ((Tj) predicts a 
mass matrix of the form (with messenger mass A^ in rji) 



kfTjerje" kfrjerje' rjTVd 
kfnWe 1 / 



vlve" 
222 



(13) 



where this matrix has to be transposed for the charged leptons. kf is the Georgi-Jarlskog 
factor which takes the values 



kf 



1 for / = d, 
-3 for / = e. 



Inserting the e suppressions of the flavon VEVs from Eq. (ITUl) the down quark mass 
matrix becomes 



e^ 



M,~ e^ e^ e'\ev,, (14) 




whilst the charged lepton mass matrix reads 

Me- I e^ -3e2 -3e^ | ewd. (15) 

£3 -3e2 

Here we assume the numerical value e ~ 0.15. Upon diagonalisation, these give mass 
ratios of e^ : e^ : 1 for the down-type quarks and y : 3e^ : 1 for the charged leptons. 
These ratios are in good agreement with quark and lepton data and also predict the 
Georgi-Jarlskog GUT scale mass relations of rrig ~ ^, f^ii ~ Stti^ and 171^ ~ rrih 
as desired. In the low quark angle approximation, the left-handed down quark mixing 
angles 9^2 ~ e, ^13 ~ e^ and 62^ ~ e^ are also predicted in agreement with data (assuming 
an approximately diagonal up quark sector which we obtain in the next subsection). 
The corresponding charged lepton mixing angles are 9I2 ~ |, Of^ ~ e^ and 62^ ~ 
3e^. Therefore, the only significant charged lepton correction to the TM mixing of the 
neutrino sector originates from 0^2 ^ -§-■, where 9c denotes the Cabibbo angle. 

3.3 Up quark sector 

Eq. ([8]) may be expanded after A^ symmetry breaking and is responsible for up-type 
quark masses 

vWe' VeVe'Ve" + VaVovlVe' vIve' \ 

vIvWe" + VaVevI'Ve' VeVo' Vd'Vl" v^ ■ (16) 

Taking the VEV hierarchy as in Eq. (TTOj) . but now adopting the messenger scale A„ ^ 
|A(i, we obtain a mass matrix with an expansion parameter e ~ 0.1, 

M„ ~ e^ e^ e^\v^. (17) 

\e' e' 1 

and an up- type quark mass hierarchy e^ : e^ : 1. This matrix gives mixing angles of 
^12 ~ ^13 ~ ^23 ~ ^^- This means that the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix is dominated by down quark mixing, except that there may be a contribution to 
^CKM fj^Qj^ f^]^Q ^p quark sector which is almost as significant as the contribution coming 
from the down-type quarks. The Cabibbo angle is still approximately 9c ~ 9f2 ~ e. 
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4 Charged lepton corrections to lepton mixing 

We have presented mixing angles which rotate the charged leptons and neutrino fields 
between the mass and fiavour bases, however these individual rotations are not what 
experiments observe. It is the combination of the two mixing matrices that appears in 
the electroweak coupling to the W boson, giving the physical mixing matrix 

UPMNS = U,,Ul^. (18) 

Here it is understood that for Lagrangians in the left right convention, Ue^ acts as 
Ue^m^m\Ul^ = {m'^^^^)'^ and U^^^ as Uy^rriyU'^^ = m!^^^. While the neutrino sector 
predicts exact TM mixing, the effect of the charged lepton corrections generates an 
experimentally detectable deviation from this in the physical parameters. In this section 
we ignore RG and CN corrections and focus only on the CL corrections. 

There are (at least) two popular ways to parameterise the PMNS matrix; firstly one 
can write t/pMNS = U23U13U12 with 

C12 Si2exp(-i(5i2) 0\ 

U12 = I -Si2exp(2(5i2) C12 , (19) 

1/ 

Ci3 Si3exp(-2(5i3)^ 

Ui3= \ 1 I , (20) 

-Si3exp{i6n) C13 

'10 

^23 = I C23 S23exp (-2(523) | • (21) 

^0 -S23exp(z(523) C23 

Here, Cij and Sij stand for cos 6ij and sinOij respectively and the 3 remaining unphysical 
phases have been rotated away, see e.g. [25]. Individual rotation matrices [4^ and U^^ 
are parameterised in the same way with relevant superscripts. The second parameteri- 
sation is that used by the PDG [16], with a Dirac phase 6 and Majorana phases a2 and 
a^; this is constructed as t/pMNS ~ R23UfP'^ R12P where the Rij are standard orthogonal 
rotations, Uf^'^ = f/13 (S13 = 6) and P = diag(l, e*~, e*~). A comparison of the two 
parameterisations, after performing a global phase redefinition to absorb remaining un- 
physical phases and obtain consistency with the convention stated in the introduction, 
shows that [25] 



S = S13- h3 - (^12, (22) 

a2= - 25i2, (23) 

as = - 2 (<5i2 + ^23) • (24) 



We can now write the parameters of f/pMNS iii terms of the neutrino mixing pa- 
rameters, with perturbative corrections from the charged lepton sector as follows [23] 
(neglecting Of^ and 6^23 as they are small) p 

S23 exp (-2(^23) ~ S23 exp (-iS^s) , (25) 

s^s exp i-^5^s) ^ 9^, exp (-15'^,) - e^.s^,, exp {-t {5!^, + S^,)) , (26) 

S12 exp {-i5i2) ^ S12 exp {-id'^^) - ^'12023012 exp {-i5l2) ■ (27) 

The dominance of the first term in Eq. fl27j) allows us to approximate 612 ~ 5]'2, while 
Eq. f l25|) gives directly ^23 ~ ^23. The phase 6is requires a more careful treatment, since 
the first term of Eq. fl2^ is larger but not dominant enough to drop the second term. 
It turns out to be possible to write 

^13 ^ ^1^3 - %^ Sin (^2^3 - 5^3 + ^12) , (28) 

^13 

assuming that -%7^ is smalllj Then the physical Dirac oscillation phase can be approx- 

"13 

imated by 

S ^ S-^^ _ ^2^3 _ 5^2 - ^ Sin (^2^3 - 5^3 + 61,) . (29) 

Turning to the resulting mixing angles, we first observe that the TM mixing of the 
neutrino sector must necessarily be a small deviation from TB mixing. We can therefore 
express our results using the neutrino TB deviation parameters [H], 



sin^r2 = ^(l + 0, sin^2^3 = i=(l + a'^), sin^r3 = ^ , (30) 

where here these parameters refer only to the neutrino sector. In terms of angles and 
phases, using Eqs. f l2^ - fl271) (see, e.g. [lU] for a discussion of this procedure), we can 
then write the TB deviation parameters for the complete lepton mixing in terms of the 
TB deviations parameters in the neutrino sector and the charged lepton corrections as, 

a ^ a", (31) 

r ^ Ir^exp (-^^^3) - 01, e^P H (<52'3 + '^12)) I , (32) 

s ^ s'' - 91, cos {6'[,- 61,). (33) 



''We note that in order to derive these equations consistently to first order, the Majorana phases 
from Eqs. (P^ -([M | must be redefined by a correction of order 9^^; this is however only a subtlety in 
the derivation and therefore we do not explicitly demonstrate this redefinition, merely point it out to 
the reader. 

^With 6*12 "^ ^, S23 "^ ^ ^iid 61^ ^ 0.15 we obtain a numerical value of -^0^ ~ |. 
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With the neutrino mixing being of TM form as given in Eq. (^, the deviation parameters 
of the neutrino sector can be shown to satisfy, see [T ^ ITT t E^. s'^ = and a^ ^ — ^ cos 6'^. 
Using this and the fact that 6*^2 ~ ^ ^-^d Eq. f l29|) . the above equations for the tri- 
bimaximal deviation parameters may be further simphfied to first order as 

a ^ — —cos 6, (34) 

r^r^-^cos(<52^3-<5r3 + 5y, (35) 

s^-ycos(5r2-5y, (36) 

again assuming that ^7^ ~ 3^ is smalL In the hmit that charged lepton corrections 
are switched off, the above results reduce to the usual TM sum rules [lll[T7l[26], s ~ 
and a ~ —^cos6. In the limit that the neutrino mixing angle 6'(^ is switched off the 
above results reduce to the usual TB sum rules [2Zj, s ~ rcosS where r ^ Oc/'i and 

The results in Eqs. flM|) - (l36|) imply the relatively simple sum rule bounds: 

\s\ < y, (37) 

|a|<i(r + ^)|cos,5|, (38) 

where, again, r, s, a are the tri-bimaximal deviation parameters, in particular r f» a/2^i3, 
6 is the CP violating oscillation phase, and 6c is the Cabibbo angle. We emphasise that 
these bounds do not include RG and CN corrections, which however are expected to 
be rather small for the case of hierarchical neutrino masses. For example, assuming 
6'i3 ~ 8° we find r ^ 0.2, and using 6'c/3 ~ 0.075 these bounds become |s| < 0.075 
and \a\ < 0.14|cos(5|. The present approximate limits from the global fit \a\ < 0.08, 
—0.06 < s < quoted in Eq. (|3]) are nicely consistent with these sum rule bounds. 

5 Conclusions 

Recent T2K, MINOS and Double CHOOZ results, together with global fits of mixing 
parameters, indicate a sizeable reactor angle ^13 ~ 8° which, if confirmed, would rule 
out TB lepton mixing. On the other hand, the TB predictions sin6'23 = 1/v^ and 
sin ^12 = l/v3 remain in agreement with global fits and continue to provide tantalising 
hints for an underlying Family Symmetry. For example, an A4 family symmetry model 
including additional fiavons in the 1' and 1" representations leads to TM neutrino mixing 
which maintains the prediction sin0i2 ~ 1/v^j at least approximately, while allowing 
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an arbitrarily large reactor angle. Indeed, as discussed in a previous paper by two of us, 
the problem in this model is in explaining why the reactor angle should be smaller than 
the atmospheric or solar angles, which follows from the fact that the additional fiavons 
would be expected to have VEVs of the same order as the other TB fiavon VEVs, with 
all undetermined coefficients being of order unity. However, apart from this drawback, 
such a model provides a simple example of a Family Symmetry model with a non-zero 
reactor angle. 

In this paper we have proposed a SUSY GUT of Flavour with a non-zero ^13 based 
on A4 Family Symmetry with additional fiavons in the 1' and 1" representations, and an 
SU{5) GUT group. The model involves an additional continuous U{1) family symmetry 
as well as three discrete symmetries designed to control the operator structure of the 
model. All fiavon representations of A4 are populated, and the main fiavon content of the 
quark sector is copied from the neutrino sector. The vacuum alignment is obtained using 
the conventional F-term mechanism. NLO terms to the mass matrices are negligible, 
demonstrating the stability of the LO matrix textures. The resulting model exhibits TM 
mixing in the neutrino sector, with the physical lepton mixing involving charged lepton 
corrections, which in turn are related to quark mixing angles. In particular, the model 
involves a Georgi-Jarlskog relation, leading to bounds on the TB deviation parameters 
\s\ < ^, |a| < ^{r + ^)|cos(5|, up to RG and CN corrections, which are in good 
agreement with current global fits. The considered model shows that it is possible to 
accommodate 613 ~ 8°, within a SUSY GUT of Favour which relates quark and lepton 
masses and mixing angles, while continuing to provide an explanation for the TB nature 
of the solar and atmospheric lepton mixing angles. 
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Appendix A: Vacuum alignment 

In order that the fiavon fields obtain the alignment presented in Eq. (^, their potential 
must be minimised in the correct way. We follow the method of [17] very closely. The 
leading order contributions to the driving superpotential are: 

Here, gi (a) = M which appears in the vacuum alignment of pJ7|; this is required since (pT 
is charged under the auxiliary symmetries and so the original structure ipj< {Mipx + ^t^t) 
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(40) 



that drives the ^px ahgnment cannot be used. Minimising with respect to ip^ gives 
The conditions from ipg and ^^ are 

(si - S2SA /sA ^ ^ /sA ^^ /S2\ M 

25-3 si - S3S1 + 5-4^ S3 + g'^u S2 + g'^u" si = , (41) 

\sl~S^S2j \S2J \sj \S,J \0j 

g^ [si + 2s2Ss) + geu^ + gju'u" = 0. (42) 

Here, (v^sj = Sj, (^) = u, (^') = u' and (^") = u". The solutions to these equations are 

g6U^ + gru'u" g'^u' + g'lu" 

u = . [AS) 




3^5 g4 

As in [19], we assume the undetermined singlets obtain their VEVs as a result of their 
soft mass parameters m^ (where s stands for singlet) being driven negative in some 
portion of parameter space. 

Appendix B: Higher order operators 

There are many higher order corrections to the mass matrices presented in Section |3] of 
this paper; most of these give negligible contributions to masses and mixings. In Tables 
E] and m we give suppressions and example terms of the NLO operators for each sector; 
it can be seen that none of these will change the LO results significantly. 
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Term 


Contributes to 


NLO Example 


c.f. LO 


FNH5 


rriD 


^^^e^e'9" ~ e6 


1 


NN 


Mr 


^le^e"fe'--e'^ 


e 


F^tH^T^ 


iMd)i2 

iMd\s 


a9' {0")^ ~ e^ 
a {9'f {e"f ~ e^ 


£6 


FipTH^T, 


iMd)u 






F^tH^T^ 


(Mdhi 

{Md)22 
{Md)23 




£4 
£3 


FipxH^T^ 


{Md)2l 

iMd)22 
{Md)23 


a^e^ {9"f ~ e9 


£3 


FiprH^T, 


iMd)si 

iMd)s2 

(Mdhs 




£7 

£6 
£ 


F^tH^T, 


(Mdh^ 

{Md),2 
iMd)s3 


a^^e' {e"f ~ e^ 


£7 
£6 
£ 



Table 3: NLO corrections in the model. The first column shows each basic term that exists 
in the neutrino, down quark (and charged lepton) Yukawa superpotential, as specified in the 
second column. A bunch of flavons is appended to these basic terms to obtain the complete 
term invariant under the symmetries. In the third column we give an example of such a bunch 
of flavons at NLO and the order of its contribution, to be compared to the LO contribution 
given in the final column. Note that in the terms contributing to M^, there is a fiavon ^pT 
already present in the basic term. It is furthermore not specified whether the LO term comes 
from an H-^ or an H-^; the reader may refer back to Eq. ([7]) if required. 
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Term 


Contributes to 


NLO Example 


c.f. LO 


TiTiH^ 


(M„)ii 


a^ee ~ e' 


e6 


TiT^H^ 


iMu)i2, iMu)21 


a^e^'^" ~ e9 


e6 


T1T3F5 


(M„)i3, (M„)3i 


a (0')^ (^")^ ~ ^ 


e3 


T2T2H^ 


{Mu)22 


{e'f 9" ~ ^ 


e3 


T2T^H^ 


(M„)23, (M„)32 


aO^e" ~ e6 


e3 


nnH^ 


(MJ33 


a^e^9'9" ~ e^ 


1 


v't 


Wo 


crVT^^^'^" ~ ?^ 


■^ 


^'s 


Wo 


a^^sOHO"fe--'^ 


^ 


e 


Wo 


aH''{9"fii"^^ 


■^ 



Table 4: NLO corrections in the model. The first column shows each basic term that exists 
in the up quark Yukawa and vacuum alignment sectors, as specified in the second column. A 
bunch of flavons is appended to these basic terms to obtain the complete term invariant under 
the symmetries. In the third column we give an example of such a bunch of flavons at NLO 
and the order of its contribution, to be compared to the LO contribution given in the final 
column. The notation 'e is simply used to denote that we are in a different sector to e or e. 

References 

[1] K. Abe et al. [T2K Collaboration], Phys. Rev. Lett. 107 (2011) 041801 
|arXiv:1106.2822j . 

[2] P. Adamson et al. [MINOS Collaboration], Phys. Rev. Lett. 107 (2011) 181802 
|arXiv:1108.0015j . 

[3] T. Schwetz, M. Tortola and J. W. F. Valle, New J. Phys. 13 (2011) 109401 
|arXiv:1108.1376j . 

[4] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, Phys. Rev. D 84 
(2011) 053007 larXiv: 1106. 6028] . 

[5] T. A. Mueller et a/., Phys. Rev. C 83 (2011) 054615 [arXiv:1101.2"663J . 



15 



[6] Y. Abe et al. [DOUBLE-CHOOZ Collaboration]. l arXiv:1112.6353 i 

[7] P. F. Harrison, D. H. Perkins and W. G. Scott, Pliys. Lett. B 530 (2002) 167 
[hep-ph/0202074] . 

[8] D. G. Lee and R. N. Mohapatra, Phys. Lett. B 329 (1994) 463 
|hep-ph/9403201|; R. N. Mohapatra, M. K. Parida and G. Rajasekaran, Phys. Rev. 
D 69 (2004) 053007 |hep-ph/03012 34]; C. Haged orn, M. Lindner and R. N. Moha- 
patra, JHEP 0606 (2006) 042 |he p-ph/0602244 |; E. Ma, Mod. Phys. Lett. A 21 
(2006) 2931 |hep-ph/0607190|; S. F. King and M. Malinsky, JHEP 0611 (2006) 
071 |hep-ph/0608021|; Y. Cai and H. B. Yu, Phys. Rev. D 74 (2006) 115005 
|hep-ph/0608022|; S. F. King and M. Malinsky, Phys. Lett. B 645 (2007) 351 
lhep-ph/0610250|; S. Morisi, M. Picariello and E. Torrente-Lujan, Phys. Rev. 
D 75 (2007) 075015 |hep-ph/0702034] ; M. -C. Chen and K. T. Mahanthappa, 
Phys. Lett. B 652 (2007) 34 [arXiv:0705.0714j; C. Luhn, S. Nasri and P. Ra- 
mond, Phys. Lett. B 652 (2007) 27 [arXiv:0706.2341]; W. Grimus and H. Kuh- 
bock, Phys. Rev. D 77 (2008) 055008 [arXiv:0710.1585j; F. Bazzocchi, S. Morisi 
and M. Picariello, Phys. Lett. B 659 (2008) 628 [arXiv: 0710.2928 j: G. Altarelh, 
F. Feruglio and C. Hagedorn, JHEP 0803 (2008) 052 [arXiv:0802.0090j: F. Baz- 
zocchi, S. Morisi, M. Picariello and E. Torrente-Lujan, J. Phys. G G 36 (2009) 
015002 [irXiv:0802.1693j; M. K. Parida, Phys. Rev. D 78 (2008) 053004 
[arXiv:0804.4571J; F. Bazzocchi, M. Frigerio and S. Morisi, Phys. Rev. D 78 (2008) 
116018 [arXiv:0809.3573j; C. Hagedorn, M. A. Schmidt and A. Y. Smirnov, Phys. 
Rev. D 79 (2009) 036002 |arXiv:0811.2"955] : H. Ishimori, Y. Shimizu and M. Tani- 
moto. Prog. Theor. Phys. 121 (2009) 769 tarXiv:0 812.5031]: F. Bazzocchi, L. Merlo 
and S. Morisi, Nucl. Phys. B 816 (2009) 204 [arXiv:0901.2086]; P. Ciafaloni, 
M. Picariello, E. Torrente-Lujan and A. Urbano, Phys. Rev. D 79 (2009) 116010 
|arXiv:0901.2236j : S. F. King and C. Luhn, Nucl. Phys. B 820 (2009) 269 
larXiv: 0905. 1686] : B. Dutta, Y. Mimura and R. N. Mohapatra, JHEP 1005 (2010) 
034 [arXiv: 091 1:22421 : S. F. King and C. Luhn, Nucl. Phys. B 832 (2010) 414 
|arXiv:0912.1344]; C. Hagedorn, S. F. King and C. Luhn, JHEP 1006 (2010) 
048 [arXiv:1003.4249]; L K. Cooper, S. F. King and C. Luhn, Phys. Lett. B 690 
(2010) 396 [arXiv:1004.3243j; H. Ishimori, K. Saga, Y. Shimizu and M. Tanimoto, 
Phys. Rev. D 81 (2010) 115009 |arXiv:1004.5"004] : S. Antusch, S. F. King and 
M. Spinrath, Phys. Rev. D 83 (2011) 013005 |arXiv:1005.0708] : G. -J. Ding, Nucl. 
Phys. B 846 (2011) 394 [arXiv:1006.4800j; K. M. Patel, Phys. Lett. B 695 (2011) 
225 [arXiv:1008.5061J; G. Altarelh and G. Blankenburg, JHEP 1103 (2011) 133 
[arXiv:1012.2697j : J. A. Escobar, Phys. Rev. D 84 (2011) 073009 [arXiv:1102.1649j . 

[9] S. F. King, Rept. Prog. Phys. 67 (2004) 107 |hep-ph/03102041; R. N. Moha- 
patra et a/., Rept. Prog. Phys. 70 (2007) 1757 |hep-ph/0510213]; R. N. Mo- 
hapatra and A. Y. Smirnov, Ann. Rev. Nucl. Part. Sci. 56 (2006) 569 

16 



[hep-ph/0603118|; C. H. Albright. [aFXiv:0905.0146l G. Altarelli and F. Feruglio, 
Rev. Mod. Phys. 82 (2010) 2701 |arXiv:1002.02TI] : H. Ishimori, T. Kobayashi, 
H. Ohki, Y. Shimizu, H. Okada and M. Tanimoto, Prog. Theor. Phys. SuppL 183 
(2010) 1 [arXiv: 1003.3552] . 

[10] S. Antusch, S. F. King and M. Mahnsky, Phys. Lett. B 671 (2009) 263 



tarXiv:071i:4727j 



|ar2£ivL0712.3759] 

larXiv:0808.2782j 



S. Antusch, S. F. King and M. Mahnsky, JHEP 0805 (2008) 066 
S. Boudjemaa and S. F. King, Phys. Rev. D 79 (2009) 033001 
S. Antusch, S. F. King and M. Mahnsky, NucL Phys. B 820 



(2009) 32 [arXiv:0810.3863j. 

[11] S. Antusch and V. Maurer, Phys. Rev. D 84 (2011) 117301 
[arXiv: 1107.3728] : D. Marzocca, S. T. Petcov, A. Ronianino and M. Spin- 
rath, JHEP 1111 (2011) 009 [arXiv:110 8.0"6T4] . 

[12] H. Georgi and C. Jarlskog, Phys. Lett. B 86 (1979) 297. 

[13] S. F. King, Phys. Lett. B 675 (2009) 347 |ar2£iv:0903.3T99]; H. -J. He, 
F. -R. Yin, Phys. Rev. D84 (2011) 033009 [arXiv:1104.2654j; Z. -z. Xing, 
Chin. Phys. C 36 (2012) 101 [arXiv:1106.3244j; N. Qin and B. Q. Ma, 
Phys. Lett. B 702 (2011) 143 [arXiv:1106.3284j; Y. -j. Zheng and B. - 
Q. Ma, Eur. Phys. J. Plus 127 (2012) 7 |arXiv:1106.4040j : S. Zhou, 
Phys. Lett. B704 (2011) 291 [arXiv:1106.4808J; T. Araki, Phys. Rev. 
D84 (2011) 037301 [arXiv:1106.5211J; N. Haba, R. Takahashi, Phys. Lett. 
B702 (2011) 388 [arXiv:1106.5926j; D. Meloni, JHEP 1110 (2011) 010 
|arXiv:1107.0"221] : S. Morisi, K. M. Patel and E. Peinado, Phys. Rev. D 84 (2011) 
053002 |arXiv:1107.069 6j: W. Chao, Y.-J. Zheng, larXiv: 1107.07381 H. Zhang, 
S. Zhou, Phys. Lett. B704 (2011) 296 [arXi v: 1107. 1097] : X. Chu, M. Dhen 
and T. Hambye, JHEP 1111 (2011) 106 [arXiv: 1107. 1589]; P. S. Bhu- 
pal Dev, R. N. Mohapatra, M. Severson, Phys. Rev. D84 (2011) 053005 
|arXiv:1107.2"378] : R. d. A. Toorop, F. Feruglio, C. Hagedorn, Phys. Lett. B703 
(2011) 447 [arXiv: 1107.3486] : W. Rodejohann, H. Zhang and S. Zhou, Nucl. Phys. 
B 855 (2012) 592 |arXiv : 1107.3970]: Q. -H. Cao, S. Khalil, E. Ma, H. Okada, 
Phys. Rev. D84 (2011) 071302 [arXiv: 1108.0570]; S. -F. Ge, D. A. Dicus and 
W. W. Repko, Phys. Rev. Lett. 108 (2012) 041801 [arXiv:1108.0964j; F. Bazzoc- 
chi, lar2£iv:1108.2497i S. Antusch, S. F. King, C. Luhn and M. Spinrath, Nucl. 
Phys. B 856 (2012) 328 |arXiv: 1108.4278] : A. Rashed and A. Datta, Phys. Rev. D 
85 (2012) 035019 [arXiv: 1109.2320]; P. O. Ludl, S. Morisi and E. Peinado, Nucl. 
Phys. B 857 (2012) 411 [arXiv: 1100393]; A. Aranda, C. BoniUa and A. D. Rojas, 
Phys. Rev. D 85 (2012) 036004 [arXiv:1110.1182]; D. Meloni, JHEP 1202 (2012) 
090 [arXiv:1110.5210j; S. Dev, S. Gupta, R. R. Gautam and L. Singh, Phys. Lett. B 
706 (2011) 168 [a rXiv:1111.1300] : A. Rashed. laDriv:1111.3072[ I. d. M. Varzielas, 

17 



JHEP 1201 (2012) 097 iarXiv:1111.3952] : R. de Adelhart Toorop, F. Feruglio and 
C. Hagedorn, Nucl. Phys. B 858 (2012) 437 (arXiv:1112.1340] : S. F. King and 
C. Luhn, arXiv:1112.1959; T. Araki and Y. F. Li, arXiv:1112.5819' S. Gupta, 
A. S. Joshipura and K. M. Patel, Phys. Rev. D 85, 031903 (R) (2012) 
|arXiv: 1 1 12.6113]: G. -J. Ding, arXiv:1201.3279; H. Ishimori and T. Kobayashi, 
lar2LiYa20 1.3429; S. Dev, R. R. Gautam and L. Singh, Phys. Lett. B 708 (2012) 
284 [arXiv: 1200753; ?• S. B. Dev, B. Dutta, R. N. Mohapatra and M. Severson, 
[^rXiv:1 202.4012l 

[14] S. F. King, Phys. Lett. B 659 (2008) 244 jarXiv:0710.0530] . 

[15] N. Haba, A. Watanabe and K. Yoshioka, Phys. Rev. Lett. 97 (2006) 041601 
hep-ph/0603116 



hep-ph/0607163 



arXiv:0809.0226 



X. G. He and A. Zee, Phys. Lett. B 645 (2007) 427 
W. Grimus and L. Lavoura, JHEP 0809 (2008) 106 
H. Ishimori, Y. Shimizu, M. Tanimoto and A. Watanabe, Phys. 
Rev. D 83 (2011) 033004 [arXw?1010.3805]; X. -G. He and A. Zee, Phys. Rev. D 
84 (2011) 053004 [ arXiv: 1 106.4359] . 

[16] K. Nakamura et al. [ Particle Data Group Collaboration ], J. Phys. G G37 (2010) 
075021. 

[17] S. F. King and C. Luhn, JHEP 1109 (2011) 042 |arXiv:1107.5"332] . 

[18] E. Ma and G. Rajasekaran, Phys. Rev. D 64 (2001) 113012 |hep-ph/0106'29l] . 



[19] O Altarelli and F. Feruglio, Nucl. Phys. B 720 (2005) 64 

[hep-ph/0504165|; G. AltareUi and F. Feruglio, Nucl. Phys. B 741 (2006) 
215 |hep-ph/0512103] . 

[20] B. Brahmachari, S. Choubey and M. Mitra, Phys. Rev. D 77 (2008) 073008 
[Erratum-ibid. D 77 (2008) 119901] [arXiv:0801.3554j; E. Ma and D. Wegman, 
Phys. Rev. Lett. 107 (2011) 061803 [ arXiv:1106.4269j . 

[21] Y. Shimizu, M. Tanimoto and A. Watanabe, Prog. Theor. Phys. 126 (2011) 81 
[ arXiv: 1105.292^ . 

[22] T. J. Burrows and S. F. King, Nucl. Phys. B 835 (2010) 174 
|ar2£ivL0909.1433j; T. J. Burrows and S. F. King, Nucl. Phys. B 842 (2011) 
107 [aFXiv:1007.2310j . 

[23] G. F. Giudice and A. Masiero, Phys. Lett. B 206 (1988) 480. 

[24] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147 (1979) 277. 



18 



[25] S. F. King, JHEP 0209 (2002) Oil |hep-ph/0204360|; S. Antusch and S. F. King, 



New J. Phys. 6 (2004) 110 |hep-ph/0405272 ; S. Antusch, S. Boudjemaa and 



S. F. King, JHEP 1009 (2010) 096 firXiv: 1003.5 498]. 
[26] S. F. King, JHEP 1101 (2011) 115 |arXiv:1011.6"T67] . 



[27] S. F. King, JHEP 0508 (2005) 105 |hep-ph/0506297|; I. Masina, Phys. Lett. B 
633 (2006) 134 [h ep-ph/0 508031] ; S. Antusch and S. F. King, Phys. Lett. B 631 
(2005) 42 [hep-ph/0508044j^ 



19 



